arXiv:1505.02879v2 [quant-ph] 11 Apr 2017 


Experimental quantum channel simulation 


He Lu/’^ Chang Liu/’^ Dong-Sheng Wang,^ Luo-Kan Chen/’^ Zheng-Da Li/’^ Xing-Can Yao,"^ Li 
Li/’^ Nai-Le Liu/’^ Cheng-Zhi Peng/’^ Barry C. Sanders/’Yu-Ao Chen/’^ and Jian-Wei Pan^’^ 

^Shanghai Branch, National Laboratory for Physical Sciences at Microscale and Department of Modern Physics, 
University of Science and Technology of China, Shanghai 201315, China 
^Synergetic Innovation Center of Quantum Information and Quantum Physics, 

University of Science and Technology of China, Hefei, Anhui 230026, China 
^Institute for Quantum Science and Technology, University of Calgary, Alberta T2N 1N4, Canada 
^Physikalisches Institut, Ruprecht-Karls-Universitdt Heidelberg, 

Im Neuenheimer Feld 226, 69120 Heidelberg, Germany 

^Program in Quantum Information Science, Canadian Institute for Advanced Research, Toronto, Ontario M5G 1Z8, Canada 

Quantum simulation is of great importance in quantum information science. Here, we report an experimental 
quantum channel simulator imbued with an algorithm for imitating the behavior of a general class of quantum 
systems. The reported quantum channel simulator consists of four single-qubit gates and one controlled-NOT 
gate. All types of quantum channels can be decomposed by the algorithm and implemented on this device. We 
deploy our system to simulate various quantum channels, such as quantum-noise channels and weak quantum 
measurement. Our results advance experimental quantum channel simulation, which is integral to the goal of 
quantum information processing. 


I. INTRODUCTION 

Quantum simulation CHS is the most promising near-term 
application of quantum computing due to the resource re¬ 
quirements for imitating some classically intractable systems 
being significantly less onerous that for other applications 
such as factorization. Experimental quantum simulation on 
closed systems is well studied using photons HI, atoms m 
and trapped ions (61. Quantum simulation of open-system 
dynamics also has variety of applications, such as dissipa¬ 
tive quantum phase transitions m and dissipative quantum- 
state engineering m, thermalization m, quantum noise gen¬ 
erators (TOl . non-Markovian dynamics idTl . and non-unitary 
quantum computing (El. Although any non-unitary quan¬ 
tum dynamics could be embedded into unitary dynamics 
over a larger Hilbert space with Hamiltonian evolution CD, 
such a direct approach entails a large computational-space 
overhead and resultant experimental complexity. Previous 
experiments have demonstrated a universal unitary-gate on 
a reprogrammable waveguide chip (T^ and have explored 
some open-system single-qubit dynamics such as quantum 
noise GsHnl , weak measurement |[I11[T3 and transpose Ei- 

Ea. 

Our approach is quite distinct from these prior achieve¬ 
ments in that our one apparatus simulates all these transfor¬ 
mations and more, in fact any single-qubit channel. We aim 
to realize a digital single-qubit channel simulator, which will 
serve as reconfigurable component of a nonunitary quantum 
circuit that would simulate nonunitary circuits. Furthermore 
our quantum simulator is “digital”. A digital quantum simu¬ 
lator is more versatile, as it is able to simulate a wider range 
of Hamiltonians O. In the sense of the analog-digital quan¬ 
tum simulation dichotomy na, for which a digital quantum 
simulator can be expressed as a concatenation of a primitive 
quantum instruction set comprising, for example, single-qubit 
gates and a controlled-not (CNOT) gate. In fact our qubit- 
channel quantum simulator requires just one CNOT gate and 


one ancillary qubit, hence is minimal in two-qubit gate cost, 
whereas ten two-qubit gates and two ancillary qubits are re¬ 
quired using standard Stinespring dilation (2^ . Our quantum 
simulator also uses far fewer single-qubit rotations than Stine¬ 
spring dilation would yield 1^ . 

In this article, we report an experimental quantum simulator 
which is imitated by a decomposition algorithm. Any single 
qubit channel can be decomposed into a mixture of two quasi¬ 
extreme channels. Experimentally, the quasi-extreme channel 
is implemented by using optical technology and the mixture 
of two quasi-extreme channels is realized by combining the 
collected data from two quasi-extreme channels. 

The article is organized as follows. Sec. [I^ provides a the¬ 
oretical introduction to channel decomposition and numerical 
simulation of the decomposition algorithm. Sec. [nl| describes 
the experimental setup. In Sec. |IV|we present the experimen¬ 
tal results, including four individual quantum noise channels 
and a simulation of weak measurement process. Finally, Sec. 
Ivl contains discussion and conclusions. 


II. THEORY 

We construct a quantum channel simulator that trans¬ 
forms a photonic qubit approximately according to any chan¬ 
nel S (211, which is a linear, trace-nonincreasing completely 
positive map that maps quantum state pto S{p). The approxi¬ 
mate experimental channel fgxp is guaranteed to map within 
a pre-specified error tolerance e with respect to o-distance 
do, which is the metric for quantifying the worst-case distin- 
guishability of the approximate from the true final state ac¬ 
cording to trace distance. Our classical algorithm for design¬ 
ing a qubit-channel simulating circuit accepts e and a 4 x 4 
matrix description of S as input and yields a description [C] 
of the photonic circuit C as output with C comprising single¬ 
qubit and a single two-qubit gate plus classical bits. A single- 
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qubit channel can be expressed as 

3 

= X! ^ij-iP'^3’ “ = (1) 

i,j=0 

with 1 the 2 X 2 identity operator and a = (X, Y, Z) the Pauli 
matrices. 

For any single-qubit channel f, p G [0,1] exists such that 

£=p£l^(l-p)£l ( 2 ) 

for each 5® a generalized extreme channel 1(241. An arbitrary 
generalized extreme channel £^ is specified by two Kraus op¬ 
erators 

Mi = Rr,{2ip)KiRrr^{26), Rr{20) := exp(-i^r • tr) (3) 
and 

f cos 3 0 \ / 0 sina\ ... 

0 c«aj’ = 0 j' 

for 0 < p < 27r. Furthermore, the Kraus operators Ki can 
be realized by the circuit shown in Fig.[^. In the circuit, 

Ry{2j) = exp(—iy 7 ) = 1 cosy — iY sin 7 , (5) 


and 


271,2 = ( 6 ) 


Each £^ has eight parameters leading to 17 parameters (in¬ 
cluding p) for arbitrary £. Random £ is generated as a two- 
qubit partial trace of a three-qubit Haar-random SU{S) ma¬ 
trix. Decomposing into Kraus operators 0 is achieved by 
guessing the 17 parameters and then optimizing by reducing 
the distance between the trial channel and the desired chan¬ 
nel £. When the trial channel £' is sufficiently close to £, the 
optimization routine terminates with the 17-parameter decom¬ 
position as output. 

We test the decomposing algorithm by numerical simula¬ 
tion. In the numerical simulation, an arbitrary channel is gen¬ 
erated from a randomly chosen unitary operator U G SU (8) 
and the channel form can be derived from Kraus operators 
Ki = {i\U\0). Five examples of input channels are 


Cl 


C 2 


C3 


C4 


C5 


/ 0.9276 -f O.OOOOi 
0.1125 -0.00391 
-0.0027-1-0.01502 

V-0.1900 +0.33832 

/ 0.7187-1-0.00002 
-0.1056 -0.10202 
0.1736 + 0.06792 

V-0.0741 + 0.26662 

/ 0.5248 + O.OOOO 2 
-0.2288 -0.15411 
-0.2444 + 0.18811 

\-0.2190-0.10141 

/ 0.3788 + O.OOOO 2 
0.1571 + 0.12111 
-0.3581 + 0.05802 

V 0.0979 - 0.0813i 

/0.4503 + O.OOOO 2 
0.1535 + 0.06041 
0.2019 -0.15411 - 

Vo. 0964 - 0 . 1834 i - 


0.1125 + 0.00391 - 
0.4846 + O.OOOO 2 
i 0.0210 -0.05541 
i 0.0787 + 0.12292 - 

-0.1056 + 0.10202 
i 0.8006 + O.OOOO 2 
-0.1944 + 0.19111 
i -0.1028 -0.14111 

-0.2288 + 0.15411 
i 0.3227 + 0.00002 
i 0.0601 - 0.0484i 
[ 0.2818 + 0.11731 

0.1571-0.12111 
0.5037 + O.OOOO 2 
i -0.2062 + 0.40691 
0.3727-0.15911 

0.1535 -0.06041 
0.4329 + O.OOOO 2 - 
-0.0169 - 0.0529i 
-0.1842 + 0.02541 - 


-0.0027 - 0.01502 - 
0.0210 + 0.05542 
0.0724 + O.OOOO 2 - 
-0.1125 + 0.00392 

2 0.1736 -0.06791 
-0.1944-0.19111 
i 0.2813 + O.OOOO 2 
i 0.1056 + 0.10202 

[ -0.2444-0.18811 
0.0601 + 0.0484i 
0.4752 + O.OOOO 2 
0.2288 + 0.15411 

-0.3581 - 0.05802 
-0.2062 - 0.4069i 
i 0.6212 + O.OOOO 2 
-0.1571-0.12111 

0.2019 + 0.15411 
-0.0169 + 0.05291 - 
0.5497 + O.OOOO 2 - 
-0.1535 -0.06041 


-0.1900 -0.3383i\ 
0.0787-0.12291 
-0.1125 -0.00391 ’ 

0.5154 + O.OOOOi / 

-0.0741 - 0.2666i\ 
i -0.1028+ 0.1411i I 
0.1056 - 0.1020i ’ 

0.1994+ 0.0000i / 

i -0.2190+ 0.1014z\ 
0.2818 - 0.1173i 
0.2288 - 0.1541i ’ 

0.6773 + O.OOOOz / 

i 0.0979 + 0.0813i \ 
i 0.3727+ 0.159H ] 
-0.1571 +0.121H ’ 

i 0.4963 + O.OOOOi / 

0.0964 + 0.1834i \ 
-0.1842 - 0.0254i ] 
-0.1535+ 0.0604i ' 

0.5671 + O.OOOOz / 


(7) 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 


Our task is to optimize the parameters specifying the chan¬ 
nel decomposition. Parametrization of unitary operator U G 
SU{2) is 

U = 


For a generalized extreme channel £{, the initial rotation 
R^{2d) is parameterized by , the final rotation 

R^{2lp) by and Kraus operators by and 

For generalized extreme channel £ 2 , the initial rotation 


cos ei-\ sin9{n ■ a) 


( 12 ) 
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is parameterized by 6 ^“^, ^, mf ^, the final rotation 

R^{2ip) by ip^^, and Kraus operators by and 

Simulation results are shown in Table 


III. EXPERIMENTAL SETUP 

Realization of Kraus operators Mi © is described by the 
circuit shown in Fig. and implemented according to the 
schematic of Fig.[^. A femtosecond pulse (150fs, 80MHz, 
780nm) is converted to ultraviolet pulses (390nm) through a 
frequency doubler LiBsOs crystal. Then the ultraviolet pulse 
(150fs, 80MHz, 390 nm) passes through two 2-mm-thick 
collinear BBO crystals, creating two pairs of photons \HV)-j 
with central wavelength of 780nm. The ultraviolet pulse 
(390nm) and the generated photons (780nm) are along the 
same direction and separated by a dichroic mirror (DM). 

The generated photons \HV)-j are separated by a PBS 
and the refiected photons are detected to guarantee that the 
transmitted photons are underway. All four photons are 
collected by the SMF and detected by the SPCM. All the 
photons are filtered by narrowband interference filters with 
AAfwhm = 2 . 8 nm (fullwidth at half-maximum) prior to de¬ 
tection. Throughout the entire experiment, the two-fold co¬ 
incidence rates for \HV )^2 are 9 x 10^s“^ and 

10 X 10^ s“^, respectively. Overall detective efficiency is ap¬ 
proximately 19%. We use a homemade Field Programmable 
Gate Array (FPGA) to record the fourfold coincidence (not 
shown here). 



FIG. 1: Experimental scheme for quantum channel simulation, a, 

The quantum circuit for Mi ([^, comprising single-qubit rotations, 
one quantum CNOT gate, and one classical CNOT gate. When the 
measurement result on ancilla is |0) (|1)), the circuit postselectively 
acts as Mo (Mi), (b) Schematic drawing of the experimental setup, 
(c) Symbols used in b. 


The single-photon rotation gates are realized by the com¬ 
bination of half-wave plates (HWPs) and quarter-wave plates 
(QWPs). The effect of HWP and QWP whose fast axes are at 
angles r and ^ with respect to the vertical axis, respectively. 


are given by the 2 x 2 matrices, 

-.j .X ( cos 2 r — sin 2 r\ 

_si„2r -cos2r)' 

C/q.p(0= 

\ —1 sm 2^ 1—1 cos 24 J 

The rotation around the y axis by angle 2ji is 

Ry{2ji) = cos 7 il -isiny^l". (14) 

The combinational operation of two HWPs set at 0° and r, 
respectively, is in the form 

C/„wp(0'>)-f/„wp(r)=f“‘f -“7- <‘5) 

\sm 2 r cos 2 r / 

We set T = 7 ^/ 2 ; then the operation is Ry{2^i). 

The rotation around 

r(T = AX^BY^CZ (16) 

by angle 20, where A, B and C satisfy + = 1, can 

be implemented by a HWP set at angle r sandwiched by two 
QWPs set at angles and ^ 2 , respectively. The combinational 
operation of the three wave plates is 

^QWp(^l)%Wp('r )^QWP(^ 2 ) 

cosO — iCsinO — sin 6 >(iA + 5)\ (17) 

— sin 0{iA — B ) cos 0 — iC AnO j 

where 

cos 0 = cos 0 cos A, A = sin 0 cos 

B = cos 0 sin A/ C = sin 0 sin A/Vn (18) 

with 

0 = ^ 1 —^ 2 , A = 2r—^ 1 —(^ 2 , ^2 = 1 —cos^ 0 cos^ A. (19) 

By appropriately choosing the angles r, and ^ 2 , rotation 
Rr{20) can be implemented. 

The two-photon controlled NOT (CNOT) gate (251 are real¬ 
ized by overlapping two photons on a polarization-dependent 
beamsplitter (PDBS). The system and ancilla photonic qubits 
are generated by shining the ultraviolet pluses on two collinear 
/3-barium borate (BBO) crystals emitting photon pairs \HV)-j 
along the pumping direction, with \H) and \ V) the horizontal- 
and vertical-polarization states, and i, j denote the path mode. 
The generated photons in the pair \HV)-j are separated by 
a polarizing beam splitter (PBS), which transmits the \H) 
component and reflects the \V) component for each pho¬ 
ton. Refiected photons 2 and 4 are collected by single-mode 
fibers (SMFs) and detected by single-photon counting mod¬ 
ules (SPCMs) to herald that photons 1 (system) and 3 (ancilla) 
are underway, respectively tZSl . 
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Cl 

C 2 

C 3 

C 4 

Cs 


Cl 

C 2 

C 3 

C 4 

Cs 

mf^ 

0.1043 

0.5254 

0.9756 

0.7444 

0.2633 


3.3938 

3.7248 

3.9264 

2.5101 

3.6886 

mfi 

0.2658 

0.7401 

0.2193 

0.1878 

0.2494 


5.3916 

2.5514 

0.7277 

1.5968 

1.6137 


0.3944 

0.1263 

0.5393 

0.9714 

0.1862 


3.1851 

6.2832 

3.2350 

4.0102 

4.0769 


0.9124 

0.9920 

0.4951 

0.2373 

0.5166 


5.6926 

4.7360 

2.3722 

1.7879 

2.5147 

mf‘^ 

0.9768 

0.0000 

0.7061 

0.8913 

0.0120 


6.1559 

3.1887 

2.2483 

0.3786 

4.5601 

^2 

^2 

0.2148 

0.9502 

0.6602 

0.0000 

0.4274 


0.1796 

3.9477 

0.5693 

1.4948 

1.1443 

nf‘^ 

0.8214 

0.2991 

0.2434 

0.7347 

0.4465 


3.4528 

1.7924 

4.7052 

4.5230 

1.7113 

nf^ 

0.5598 

0.4585 

0.7350 

0.4100 

0.8948 

P 

0.4675 

0.7001 

0.8370 

0.9032 

0.4922 


3.8118 

4.7982 

6.1211 

5.8303 

5.4589 

e 

0.0031 

0.0009 

0.0040 

0.0006 

0.0045 


TABLE I: Numerical simulation of the channel decomposition for five randomly chosen input channels and channel Ct- Error e is the actual 
error from the simulation. 
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XI 

XX 

XY 

XZ 
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YY 

YZ 

Z I 

ZX 

ZY 

ZZ 


,-X>-N _X>-N_X>^N 
■ XXXX XXXXNNNN 


U.zt) 

I 

I 


.xxn-x^'^-x^n_xxn 
--XXXXXXXXNNNN -0-25 


EIG. 2: Process matrix of the CNOT operation, a. Real elements of ideal process matrix (imaginary elements are identically zero), b. The 
real elements of Xexp. c. Imaginary elements of Xexp. 


We experimentally characterize the quantum CNOT gate 
via quantum process tomography (QPT) technology ll27l and 
obtain gate fidelity F = 0.83 ± 0.02 as shown in Fig|^ 

The classical CNOT operation is a classical logic operation 
that flips the system-qubit state 1' conditioned on the mea¬ 
surement result of ancillary qubit 3'. Experimentally, clas¬ 
sical CNOT is effectively statistically simulated: we set the 
measurement basis of ancillary photon 3' on \H) or \ V) with 
equal probability. No further operation on system qubit 1' oc¬ 
curs when the measurement basis choice of ancillary photon 
3' is \H), whereas an X operation (an HWP set at 45°) is 
applied on the system qubit 1' when the measurement basis 
choice of ancillary photon 3' is \ V). If the ancilla-qubit mea¬ 
surement result is |iT) (|y)), the simulator is described by Mq 
(Ml). 

The probability p is also statistically simulated. We first set 
up the circuit for simulating and collect data for time ti. 
Then we convert the circuit to the case of simulating 
collect data for time t 2 - Combining these data yields E{p) = 
{tiSi + ^ 2 ^ 2 )/(^i + ^ 2 )- By choosing ti and ^2 appropriately, 
any p e [0,1] can thus be simulated. 

We emphasise here that the rotations in our experiment 
are realized manually, and the classical CNOT gate is imple¬ 
mented by inserting an HWP according to the projector choice 
on photon 3'. The mixture of f ° and £2 is statistically simu¬ 


lated by collecting data from and £2 with different times. 
In our experiment, simulating one single-qubit channel needs 
us to run the setup four times and then combine the collected 
data. 


In fact, the four runs can be embedded into one run. Here, 
we also propose an experimental scheme to simulate any 
single-qubit channel within one run. Below we summarize the 
experimental scheme. As shown in Fig.[^ a beam of light with 
two different colors, with central wavelength of Ai and A 2 re¬ 
spectively, is split by a dichroic mirror(DM). The transmitted 
color Ai go through f ° as shown in Fig.[^ The reflected color 
Ai first go to an automatically controlled attenuator and then 
go through f The attenuator setting depends on the parame¬ 
ter p G [0,1]. For example, the attenuator block half of A 2 is 
equivalence of simulating p = 2/3. Finally, Ai and A 2 is re¬ 
combined on another DM and injected to the detector. Ah the 
rotations can be realized by an automatically controlled wave- 
plate and the classical CNOT can be replaced by feedforward 
technology. 
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DM 


s: 


DM 




attenuator 


IV. RESULTS 


A. Random channel 


FIG. 3: Proposed experimental scheme to simulate single-qubit We first show that a randomly chosen channel is accurately 
channel with one run. simulated with our setup. Randomly chosen input channel £ 


/ 0.3938 

0.0075 - 0.0739i 
0.0172 - 0.0155i 
\0.0801 + 0.0614i 


0.0075 + 0.0739i 
0.1594 

-0.0733 + 0.0801i 
-0.066 - 0.0172i 


0.0172 + 0.01551 
-0.0733 - 0.0801i 
0.2241 

-0.014 + 0.0751 


0.0801 - O.O 6 I 4 A 
-0.066+ 0.0172i 
-0.014 - 0.075i 
0.2228 ) 


( 20 ) 


mf^ 

0.1896 


0.7948 

nf^ 

-0.7813 

nfi 

0.5804 

mf" 

0.3901 

mp 

0.5051 


-0.0919 


0.9817 


0.187r 

p£i 

0.267r 


0.847r 

^£2 

0.407r 


0.427r 


0.367r 


-0.757r 


0.567r 

p 

0.6 


TABLE II: Decomposition of channel £ 


is realized by the circuit in Fig.[^ with appropriate parameters 
specified. The decomposition of £ is shown in Table |I^ For 
p = 0.6, we set + = 60s and ^2 = 40s. 

To verify accurate channel simulation, we use QPT to re¬ 
construct the matrix representation of £. Figure shows the 
experimentally reconstructed fexp matrix. We calculate the 
process fidelity 

Fp = IV (21) 

between the reconstructed matrix fgxp and £ and discover 
Fp = 0.94 ± 0.02. 



FIG. 4: Reconstructed process matrix of the randomly chosen chan¬ 
nel in Eq. ( [^ . (a) The real part of the ideal channel £. (b) The 
imaginary part of the ideal channel £. (c) The real part of the ex¬ 
perimentally constructed channel ^exp. (d) The imaginary part of the 
experimentally constructed channel ^exp. 


Average fidelity is 1^ 

F = (2Fp + l)/3 = 0.96 ± 0.01. (22) 

As further analysis, we calculate the trace distance 
D{£, fexp) = Ti\£- fexp|/2 = 0.22 ± 0.02. 


( 23 ) 
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Fidelity Fp is related to D by the inequality lua 

1 - < D < x/l - ^’p- (24) 

In our case (Fp = 0.94), the upper and lower bounds of D are 
0.06 and 0.24. 

B. Amplitude Damping Channel 


A 

£t 

P 

a 

p 

7i 

72 

7?m(2(5) 

Rn{2ip) 

0 

0 

0 

7r/4 

-7r/4 

none 

none 

1 

0.36 

O.IOStt 

O.IOStt 

7r/4 

-0.04471 

none 

none 

1 

0.5 

7r/4 

7r/4 

7r/4 

0 

none 

none 

1 

0.75 

7t/3 

7r/3 

7r/4 

7r/12 

none 

none 

1 

1 

7tI2 

7r/2 

7r/4 

7r/4 

none 

none 

1 


The amplitude damping (AD), or decay channel can be de¬ 
termined by two Kraus operators 


Ko = 



,Ki = 


Cf) 


(25) 


Table |n^ shows the setting of parameters of the 
simulator. Fig. |^(e) shows the geometric interpre- 


TABLEIV: Parameters of bit-flip channel for the given A. 


For different A, the setting of the parameters are shown 
in Table [V] Fig. [^(g) shows the geometric interpreta¬ 
tion of ideal(experimental) phase-flip channel for A G 
{0,0.36,0.5,0.75,1}. 


A 

£l 

P 

a 

p 

7i 

72 

7?m(2<5) 

Fn(2p) 

0 

0 

0 

7i/4 

-7i/4 

none 

none 

1 

0.36 

0.10371 

0 

0.271 

-0.271 

none 

none 

1 

0.5 

7i/4 

0 

371/8 

-371/8 

none 

none 

1 

0.75 

7i/3 

0 

7i/12 

-7i/12 

none 

none 

1 

1 

7i/2 

0 

0 

0 

none 

none 

1 


TABLE III: Parameters of amplitude damping channel for the 
given A. 


E. Depolarizing channel 

The depolarizing channel, which is known as a white-noise 
channel, has the form 

£{p) = {l-\)p+^{XpX + YpY + ZpZ). (28) 

o 

The setting of the parameters are shown in Table[Vl| Fig.|^(h) 
shows the geometric interpretation of ideal(experimental) de¬ 
polarizing channel for A G {0,0.36,0.5,0.75,1}. 


tation of ideal(experimental) AD channel for A G 

{0, 0.36, 0.5, 0.75,1}. F. Weak measurement 


C. Bit-flip channel 


The bit-flip channel has two Kraus operators in the form 


Ko = VT^l = ^ 


0 




Ki = VXX = 


0 

,^/A 0 ; 


(26) 


For different A, the corresponding channel parameters are 
shown in Table Fig. [^(f) shows the geometric inter¬ 
pretation of ideal(experimental) bit-flip channel for A G 
{0,0.36,0.5,0.75,1}. 


D. Phase-flip channel 

The phase-flip channel has two Kraus operators in the form. 


Ko = Vl - A1 = 

Ki = \/A2’ = 




0 


0 ’ 

Vx 0 \ 

0 -Vx ■ 


(27) 


We show that our apparatus successfully simulates trace- 
decreasing channels, such as weak measurement followed by 
measurement reversal, which is a strategy for offsetting ampli¬ 
tude damping Sad (25) at the cost of losing particles through 
postselection csiiiMa. For a single-qubit input state p, this 
strategy is 


p^ p' = N[£ad{MpM^)]N^ (29) 


for weak measurement M = diag(l, v^l — pi) with pi ^ 0 
and weak measurement reversal N = diag(v^l — p 2 , 1) with 
P 2 ~ 0. A successful outcome corresponds to high fidelity 

/ ,_X 2 


Tr 


with success probability Tr p'. 


Larger pi corresponds to superior protection and smaller 
success probability. Seeking to explore the trade-off between 
success probability and weak measurement strength pi, we 
choose a fairly strong measurement strength pi = 0.8 and 
then let p 2 = Pi + A(1 — pi) if damping parameter A is 
given (this relation between pi and p 2 is the “optimal strat¬ 
egy”); otherwise p 2 = pi if A is unknown (the “non-optimal 
strategy”) Ool. Theoretical and experimental state fidelity 
results for input states |F), |F), \D) = l/x/2{\H) + \V)) 
and \L) = l/x/2{\H) + i |V)) are shown in Fig.for three 
cases: pure amplitude damping, non-optimal measurement 
strategy, and optimal measurement strategy. Note that \H) 
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(e) (f) (g) (h) 

FIG. 5: The geometric interpretation of the quantum-noise channels on a Bloch sphere, (a)-(d), the ideal geometric interpretations of amplitude 
damping channel, bit-flip channel, phase-flip channel and depolarizing channel with A = 0(blue dots), A = 0.36(red dots), A = 0.5(green 
dots), A = 0.75(purple dots) and A = 1 (yellow dots), (e)-(f), the geometric interpretations of experimentally reconstructed channels. 


A 



P 

a 

/? 

71 

72 

Rn,{5) 

RnW 

a 

13 

71 

72 

Rn,{S) 

Rr,W 


0 

0 

0 

7r/4 

- 7 v /4 

none 

none 

TT 

0 

- 7 t /4 

7 v /4 

none 

none 

1 

0.36 

0 

0 

7 v /4 

- 7 v /4 

none 

none 

TT 

0 

- 7 v /4 

7 v /4 

none 

none 

0.64 

0.5 

0 

0 

7 v /4 

- 7 v /4 

none 

none 

TT 

0 

- 7 v /4 

7 v /4 

none 

none 

0.5 

0.75 

0 

0 

7 v /4 

- 7 v /4 

none 

none 

TT 

0 

- 7 v /4 

7 v /4 

none 

none 

0.25 

1 

0 

0 

7 v /4 

- 7 v /4 

none 

none 

TT 

0 

- 7 v /4 

7 v /4 

none 

none 

1 


TABLE V: Parameters of phase-flip channel for the given A. 


input is immune to Sad, but, due to experimental imperfec¬ 
tion, the fidelity for input state \H) is not exactly 1. We find 
that the optimal strategy provides the best protection, and the 
experimental results agree with the theory for all three cases. 


V. CONCLUSION 

In this article, we demonstrate that a digital channel sim¬ 
ulator can be realized via linear optics. Any open-system 
quantum dynamics and quantum channels on single qubit can 
be simulated in our system. For multi-qubit channel simu¬ 
lation, decomposition algorithm has been extended to qudit 
channels ISTl . In large-scale channel simulation, linear op¬ 
tics system might be retarded by the probabilistic CNOT gate. 
However, other systems, such as superconducting qubit and 
trapped ions, can benefit from our results. Our demonstration 
can serve as a foundation for future experimental simulations 


employing networks of qubit channel simulators. Such net¬ 
works could serve to simulate general dissipative many-body 
dynamics including the interplay between dissipative and uni¬ 
tary processes lf32l and dissipative universal quantum compu¬ 
tation |[T^ by combining two-qubit entangling gates with the 
qubit-channel simulators. 
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a 
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p 

71 
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0 
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0 

Y 
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7r/4 

7r/4 
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0 
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7r/4 
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0 
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7r/2 
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7r/4 

0 

Y 
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TABLE VI: Parameters of depolarizing channel for the given A. 
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(c) (d) 



FIG. 6: State-fidelity dynamics for the amplitude-damping channel with weak measurement protection. Dashed lines represent theoretical 
results, and dots represent experimental results. Red, blue and black colors represent no weak measurement, non-optimal measurement 
strategy (pi = P 2 ) and optimal measurement strategy (p 2 = pi + A(1 — pi)), respectively. The input state is (a) \H), (b) \V), (c) \D), and 
(d) \L) 













